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On the linear level elliptic equilibria of Hamiltonian systems are mere superpositions of harmonic oscilla- 
tors. Non-linear terms can produce instability, if the ratio of frequencies is rational and the Hamiltonian 
is indefinite. In this paper we study the frequency ratio ±1/2 and its unfolding. In particular we show 
QQ that for the indefinite case (1:— 2) the frequency ratio map in a neighbourhood of the origin has a critical 

r- H point, i.e. the twist condition is violated for one torus on every energy surface near the energy of the 

equilibrium. In contrast, we show that the frequency map itself is non-degenerate (i.e. the Kolmogorov 
| | ' non-degeneracy condition holds) for every torus in a neighbourhood of the equilibrium point. As a by 

{fy product of our analysis of the fequency map we obtain another proof of fractional monodromy in the 

1: — 2 resonance. 

•i-H 

l^j 1 Introduction 

The problem of resonant equilibria in Hamiltonian systems is quite intricate. The set of resonant frequencies 
is countable and dense within all frequencies. It is therefore generic for a Hamiltonian system to have all 
equilibria non-resonant, but the set of all such Hamiltonians is not open. An approach that allows to recover 
an open and dense set of Hamiltonians is to make a distinction between resonances of high and low order. The 
latter concern equilibria with local dynamics significantly different from that of any Birkhoff normal form. 
In two degrees of freedom this makes the 0:0, 0:1, 1:±1, 1:±2 and 1:±3 resonances low order resonances. For 
high order resonances the local dynamics is essentially governed by a Birkhoff normal form, cf . [3D] . 
(~ — Thus it requires an external parameter A to encounter a resonant equilibrium in a Hamiltonian system. 

The role of A is to detune the frequencies. In fact, the 0:0 resonance — an equilibrium with nilpotent 
linearization — is of co-dimension 2 and needs generically one detuning parameter for each of the two 
frequencies, cf. [5]. The 1:1 resonance is special as well : already the linear centraliser unfolding needs 
3 detuning parameters. In [7 a universal unfolding of the full (nonlinear) problem with 7 parameters 
is computed. The 0:1 and 1:— 1 resonances lead under parameter variation to hyperbolic eigenvalues. This 
triggers the normally elliptic centre-saddle (or Hamiltonian saddle-node) bifurcation [4] and the Hamiltonian 
Hopf bifurcation |33j , respectively. 

In this paper we consider both the 1:2 and the 1:— 2 resonance ; we use the notation a — ±1 in order to 
treat these cases in a parallel fashion. Moving the resonant equilibrium to both the origin (q,p) — of the 
phase space M 4 and the origin of parameter space to A = 0, the linear part of the Hamiltonian vanishes and 
the quadratic part at the resonance reads 

H 2 (q,p) = | [p\ + q\ + 2a{p 2 2 + ql)] . 
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The eigenvalues are ±iw, ±2iw and thus in 1:±2 resonance. We scale the frequency uj to 1 by a linear change 
of time. The linear centraliser unfolding 

G 2 (g,p;A) = l -{ P \ + ql) + a( P 2 2 + q 2 2 ) + ^{pI + qI) 

consists in adding a convenient detuning term. 

Our interest lies in the full nonlinear dynamics defined by a family H(q,p;X) of Hamiltonians with 
quadratic part G 2 . In addition to the linear terms we also set the constant term equal to zero. Thus we are 
confronted with 

H(q,p;X) = G 2 (q,p;X) + G 3 (q lP ;X) + G 4 (q,p;X) + ... (1.1) 

where Gk is a homogeneous polynomial of order fc in (q,p) that may have A-dependent coefficients. Close 
to the origin Gk dominates Gi for k < I and we claim that already the cubic polynomial 

H 3 (q,p;X) = G 2 (q,p;X) + G 3 (q,p;X) 

contains the desired information, provided that the average G 3 of G 3 over the periodic flow of Xjj 2 does not 
vanish at A = 0. Indeed, our first step will be to derive a normal form 

H 3 (q,p;X) = G 2 (q, P ;X) + G 3 (q, P ;X). (1.2) 

Note that this would also be the first step to compute a Birkhoff normal form, for which the terms of odd 
order vanish. Because of the resonance we have an unremoveable odd order term. The result is called 
resonant Birkhoff normal form. 

In case G 3 does vanish identically we may proceed to compute the Birkhoff normal form up to order 
4, and for non-degenerate quartic terms this special 1:±2 resonance becomes a high order resonance. The 
tricky bit is when G 3 does vanish at A — 0, but not identically in A. This case is of co-dimension 2 and one 
would then use a coefficient c in G 3 as second unfolding parameter. 



By construction the normal form Hamiltonian ( 1.2 ) is invariant under the periodic flow generated by H 2 . 



Our second step consists in the reduction of this ^-symmetry. The origin (q,p) — is fixed whence S 1 does 
not act freely and we have to resort to singular reduction [9J. This reduction is performed by choosing a 
Hilbert basis {tt\, tt 2 , ir 3 , ka}, i-s. generators of the ring of S* 1 -invariant functions as (global) co-ordinates. A 
certain drawback is that this basis is not free, whence we have to work on the 3-dimensional semi-algebraic 
subvariety of ]R 4 that is defined by 

^ = n 2 3 + ill, 7T 1 >0,7r 2 >0. (1.3) 

The normal form H 3 induces a Hamiltonian TL on this reduced phase space. Being of one degree of freedom the 
corresponding reduced Hamiltonian system is readily analysed. We do this both globally on the symplectic 
leaves of the reduced phase space and locally in well-chosen co-ordinates. See Elipe et al [TSj HH] for a 
similar approach in the setting of stability analysis. Due to the geometry of the reduced phase space, it is 
convenient to plot the energy surface in the space of the invariants (figures [2] and [3]) . This technique was, in 
fact, already known to Alfriend [1] and even earlier to Breakwell and Pringle [2]. 

Our third step concerns the dynamics in two degrees of freedom. This consists of two parts. First we 
reconstruct the dynamics defined by the normal form H 3 from the reduced Hamiltonian TL. In the case 
a = +1 of the 1:2 resonance this makes K 4 a ramified torus bundle, with regular fibres T 2 reconstructed 
from the periodic orbits of the reduced system. In the case a = — 1 of the 1:— 2 resonance the reduced system 
always has unbounded motion. It may also have bounded orbits depending on the detuning. We show by a 
direct computation that R 4 becomes a ramified fibre bundle, where the regular fibres are 2-tori and cylinders 
reconstructed from bounded and unbounded orbits, respectively, in the pre-image of the regular values of 
the energy-momentum mapping. In order to give a global overview of the dynamics we then go on and 
construct in both cases a = ±1 the sets of critical values of the energy-momentum mappings. 



The fourth step relates the reconstructed integrable dynamics defined by (1.2) to the "original" dynamics 



defined by (1.1 1. For a = +1 this perturbation analysis is straightforward. Things turn out to be more 



interesting for the 1:— 2 resonance. Indeed, we show that the ramified fibre bundle defined by (1.2 1 displays 
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fractional monodromy around the singular fibre containing the origin (q,p) — (in the case a = +1 the 
corresponding singular fibre consists of the origin alone). This phenomenon is treated in Section [5] The 
main technical issue is to define the notion of frequency and frequency ratio (which shows monodromy) for 
the non-compact cylinders. The approach taken in |T7] is to add higher order terms to the Hamiltonian 
that make the energy surface compact. We follow a different approach along the lines of [25] where only 
the singular part of the dynamics near the equilibrium point is considered. This singular part captures the 
dynamics up to a smooth (but unknown) function that encode the influence of higher order terms. However, 
the singular terms dominate the asymptotic behaviour of the frequency map when sufficiently close to the 
equilibrium, and thus nothing needs to be known about the higher order terms, except their existence and 
smoothness. 

The unfolding shows that our situation is different from the unfolded 1:— 1 resonance. When passing 
through a 1:— 1 resonance, (non-fractional) monodromy is created in a supercritical Hamiltonian Hopf bi- 
furcation and local monodromy is turned into non-local [34] island |16| monodromy in the subcritical case. 
The situation is different for the 1:— 2 resonance (where no distinction super/subcritical exists). Fractional 
monodromy |28] is local at A = and turns into island monodromy for both positive and negative values 
of A. Stated differently, the system has fractional monodromy before and after resonance. When passing 
through a 1:— 2 resonance island monodromy momentarily contracts to local monodromy. 

In the final section we come to our main new results on the 1:— 2 resonance. We show that the isoenergetic 
non-degeneracy (or twist) condition is violated for a single torus on every energy near the critical energy. 
This shows that the principle "low order resonance leads to vanishing twist" (already established for low 
order resonant fixed points of maps [H1H2]) ' s a l so valid for low order resonant equilibria. This shows that 
the principle "monodromy causes vanishing twist" , which has been established for the Hamiltonian Hopf 
bifurcation (13] and for focus-focus equilibria in general |15] also holds in the more degenerate (and more 
difficult) situation of the 1:— 2 resonance. The second main new result is that the frequency map for tori 
near the 1:— 2 equilibrium is non-degenerate, i.e. the Kolmogorov condition for the KAM theorem holds. 
Again, this is similar to the situation for focus- focus points [TS]. The main difference with the focus-focus 
equilibrium is that this equilibrium was assumed to be non-degenerate (in the sense that the linearisation of 
the commuting vector fields generate the Cartan sub-algebra); whereas in the 1: — 2 equilibrium it is always 
degenerate. Our analysis shows that this implies that the period only has algebraic blowup, as opposed to 
the usual logarithmic one. 

Our results are valid for all completely integrable systems that agree up to the cubic terms with the 
normal form we analyse. The higher order terms do not change the asymptotic properties of the frequency 
map when sufficiently close to the origin. An arbitrary Hamiltonian system near the 1 : —2 resonance is 
generically non-integrable due to higher order terms. Thus is does not posses a frequency map. Using 
our results about the non-degeneracy of the frequency map, then KAM theory can be invoked to show the 
persistence of tori whose frequencies are Diophantine. Moreover, monodromy can be defined for the nearby 
non-integrable system [3J. 

2 Normal form 

The art of normalizing Hamiltonian Systems has been periodically (re)discovered and forgotten and goes at 
least back to Birkhoff. For treatments similar to ours see e.g. [33 ] 15] [13 ] . The resonant quadratic part Hi 
generates the S^-action 



q : 



^xC 2 

(t,Zl,Z 2 ) 



C 2 

(zie 1 *, z 2 e 2<rit ) 



where Zk 



(pk + igfe)/v2, k = 1,2. Introducing 




this action has the invariants 



7Ti = zizi , 7r 2 = z 2 z 2 , tt 3 = Re(zfz% ) , 7T 4 = Im(z 2 ^2 ) . 
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These form a Hilbert basis of the ring (C(R 4 )) e of g-invariant functions, i.e. given a g>-invariant function 
/ : M 4 — ► R there exists g : M 4 — ► K such that 

/(<?,?) = 9{n{q,p)) ■ 

Normalization of H(q,p; A) with respect to this action gives a Hamiltonian that is a function of the invariants 
only. The only cubic polynomial that is a polynomial of the invariants is the linear combination air 3 + 6^4. 
A simultaneous rotation in both the (gi,pi)-plane and the ((72,^2) -plane turns this expression into C7r 3 with 
c 2 = a 2 + b 2 . Unless c = we assume that c > 0. This yields the following result [TT| [25j. 



Theorem 1. Let the family H(q,p; A) of two-degree-of-freedom Hamiltonian systems be given by (1.1 1, i.e. 
with quadratic part unfolding the 1:±2 resonance. Then there is a symplectic co-ordinate transformation %j) 
that turns H(q,p; A) into normal form 

(H o tp) (q,p; A) = G2(q,p;X) + G 3 (q,p;X) + higher order terms 

with G 3 (q,p; A) = c(A) [{p\ - q\)p 2 + 2<rp 1 q 1 q 2 \ . 

Introducing the scaling (g, p, H) 1— > (eg, ep, e 2 _ff ) zooms into the origin and reveals the smallness of the higher 
order terms as the normal form reads 

(Hoyj)(q,p;X) = G 2 (q, P ;X) + eG 3 (q,p;X) + 0(e 2 ) . 

Restricting to a small interval ]— Ao, Ao[ in parameter space defined by 

c(A) = c(0) + 0{e) 

we may furthermore put c(A) 



c(0) 



perturbation of (1.1 1. Note that if c'(0) 7^ 0, then Ao 
This turns ( |l~2| into 



: c in G3 whence the normal form (1.2 1 turns out to be an e 2 -small 

eA as well. 



O(e), Therefore we may rescale A 



H 3 (q,P] A) 



I ^ (Pi + Qi) + c [(Pi ~ ii)P2 + 2o-pigig 2 ] 



From now on we assume that the non-degeneracy condition c > holds. The additional scaling (A, e) > 
(Ac, c _1 e) would yield c = 1. We keep c in the calculation so that we can absorb additional constants when 
they appear. Moreover, later on we want to (re)-introduce the scaling in a slightly different way in the image 
of the energy-momentum map. 



3 Reduction 

The normalized truncated cubic Hamiltonian is Liouvillc integrable with integral iJ 2 - The flow generated by 
H2 is periodic, and hence it is the momentum of an ^-action. Singular reduction of this symmetry gives a 
system with one degree of freedom. Instead of working with tti , ~ki , ix 3 , 7r4 it is convenient to use a slightly 
different Hilbert basis of the ring (C(K 4 )) e of g-invariant functions that contains the generator of the S 1 
action. Therefore we drop 7T2 in favour of 

r, = tti + 2a7r 2 = \{pl + ql) + a{p 2 2 + ql) 

Whence the relation between the invariants reads R^^tti, ir 3 , n^) — with 

i?,,(7ri,7T3,7r 4 ) = ^7rf(77 — tti) — (7r| + 7r|); 

while the inequalities tti > 0, 7T2 > turn into an interval 2^ of admissible values for ~k\. In the case a = +1 
of the 1:2 resonance this interval is bounded and given by 
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Figure 1: The reduced phase spaces that occur in the 1 : ±2 resonance. Upper left: a = +1 and h% > 
(not shown: = for a = +1 and fi2 < and P^ 1 " = {0}); upper right: a = — 1 and > 0, lower left: 
a = — 1 and /12 < and lower right: cr = — 1 and /12 = 0. 
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Figure 2: Section 7r 4 = and 7T1 + 27r2 = /12 of the reduced phase space for a = +1, h% = 1, together with a line of 
constant Hamiltonian h = hi + A71-1 + C7T3 = 1.2, A = 1/2, c = 1. 

whereas for u = — 1, the invariant 77 may assume negative values as well and 
I~ — [max{0, 77} , 00 [ 

is unbounded. The Poisson bracket relations between the invariants are easily computed and given in tablc[T| 
where 

/On) = |(2bm-3*?) = g?a. 

Because 77 coincides with Hi it was to be expected that this invariant is a Casimir of the reduced Poisson 
structure. 



Table 1: The reduced Poisson bracket. 
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Fixing the value ft.2 of 77 the one-degree-of-freedom dynamics takes place in the surface of revolution (see 
figure [lj 

P£ = { (7ri,7T 3 ,7r 4 ) e M 3 I i? h2 (7rx,7r3,7r 4 ) = 0, 7Ti € } 
with Poisson structure 

{f, 9} = (V/ x V. 9 I VJfc) . 

Note that is compact and is unbounded. The conical singular point at 7n = (unless a — —1 and 
/12 > 0) comes from the isotropy subgroup Z2 of the 5' 1 -action p at (0, Z2). When cr = — 1, such points satisfy 
/z 2 < 0. 



G 



^3 




Figure 3: Section m = and m — 2n2 = of the non-compact reduced phase space for a = — 1. For h% = — 1 
(top), there is a singular point at the origin. The line of constant Hamiltonian is determined by Ah — 0.2, A = 1/2, 
c = 1. For = 0.1 (bottom), the reduced phase space is smooth. The line of constant Hamiltonian is given by 
Ah = 0.0565, A = 1/2, c = 1. 
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The reduced Hamiltonian is obtained by expressing if 3 in the basic invariants and reads 

H 3 (7ri,7r3,7r 4 ; A) = 77 + e {Xiri + cir 3 } , (3.4) 

where the factor 2V2 has been absorbed in the factor c. Subtracting the constant value hi of 77 and rescaling 
time by e yields 

H(7ri,7r 3 ,7r 4 ; A) = Xtti + C7r 3 . 

It is important to keep in mind that subtracting hi is only a valid operation in the reduced system, but not 
in the full system. Altering the Hamiltonian by a constant in general makes no difference to the dynamics. 
Eventhough H 2 is constant along trajectories of the full averaged system, it does contribute leading order 
terms to the equations of motion, and thus influences the frequencies and rotation number. In the reduced 
system the function 77 (with value hi) is a Casimir of the reduced bracket. Hence removing it does not change 
the equations of motion. So even though in the unreduced system we should not remove Hi, we will see 
that the fractional monodromy is created by the non-linear terms G 3 , and the monodromy is not changed 
by removing Hi. 

The reduced integral curves are given by the intersections 

where Ah = h — /12 is the difference between the value h of the energy and the value hi of 77, i.e. the value 
of the cubic terms in the Hamiltonian. In addition to the singular points, we have regular equilibria where 
Ti.~ 1 (Ah) is tangent to P% . Because the latter is a surface of revolution and the former is a plane parallel 
to the 7T4-axis this requires tt^ = 0. Hence, we merely have to compute the points where the straight line 

Ah A 

7T 3 = - -7Tl 

C C 

is tangent to the cubic curve when 7Ti S Z£ and 

13 = |fl"i( } 7- 7r i)> ( 3 - 5 ) 

see figures [2] and [3j When this happens at the singular point iri — 0, a Hamiltonian flip bifurcation occurs 
as hi is varied, cf. |21j . In the full system this corresponds to a supercritical period-doubling bifurcation. 
This explains the phase portraits when a = +1. For a = — 1 the subcritical Hamiltonian flip bifurcation at 
the singular equilibrium for hi < is accompanied by a centre-saddle bifurcation (of regular equilibria) for 
h 2 > 0, both emanating from the non-conical singularity tt\ — at A = 0. This is best illustrated in terms 
of the critical values in the image of the energy-momentum map, see below. 

3.1 Rational Parametrisation 

An interesting feature of the 1:±2 resonances is that the two semi-algebraic varieties and PZ h together 
form the cubic R^(0). Only when h 2 < the reduced phase space P^ 2 is empty and the cubic (0) has 
an additional isolated point at the origin. Since this cubic always has a singular point there exists a rational 
parametrisation given by 

(7Ti,7r 3 )(s) = (h 2 - 2crs 2 ,S7Ti) 

where s = 7r 3 /7ri is the slope of the line through the origin that intersects the cubic at (7Ti, 7r 3 ). Similarly there 
is an explicit polynomial parametrisation using the inhomogeneous co-ordinates Si = 7r 3 /7ri and si = ^/Vi 
such that 

(7ri,7r 2 ,7r 3 ,7r 4 )(si,S2) = (h 2 - 2a(sf + s\), s\ + s|, si7Ti, s 2 iti) . 

This stereographic projection from the singular point onto the plane ni = 1 gives co-ordinates on the 
reduced phase space. Here (711,^3,^4) are like homogeneous co-ordinates, while si,si are inhomogeneous 
co-ordinates on the plane ~k\ — 1. A similar construction can be done with the plane ni = \. 
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-1 -0.5 0.5 1 

Si 

Figure 4: Phase portrait of H(s; A) for a = ±1, h 2 = ±1, A = 1/2, c = 1. 



Theorem 2. The inhomogeneous co-ordinates (sj., S2) = (^3/^1, ^4/^1) are symplectic co-ordinates (up to 
a multiplier) on the singular reduced phase space: 

{si,s 2 } = | . 

The reduced Hamiltonian TI(tt; A) = Xiti + CK3 in these co-ordinates is 

H(s;X) = (X + csi){h 2 - 2a{s\+sl)) . 

Proof. By direct calculation. □ 

The level lines of 7i(s; A) define the planar phase portrait, see figure [4] for an example. The singular point 
of the reduced phase space (in 7Tj) is blown up to the invariant circle h 2 — o~(s\ + s 2 ) in the phase portrait. 
If the line si = — A/c intersects the circle it contains two equilibrium points. For a = 1 the relevant part 
of the phase portrait is the closed disk inside this circle, while for a = — 1 and h 2 < it is the closure 
of the complement. Recall that the two reduced phase spaces for ah 2 > are just different parts of the 
same singular cubic curve, and in the same way the union of the two phase portraits gives all of the plane 
M 2 = { Sl ,s 2 }. 

When there is no singular point in the reduced phase space, ah 2 < 0, the phase portrait is all of 
R 2 = {si, s 2 }. The circle seen in the case ah 2 > now has negative radius. 

The other prominent feature of the phase portrait is the invariant line si = —A/c. It is the image of the 
energy surface Ah = 0, which is a plane through the origin = \tt\ + c7T3 intersected with the reduced phase 
space. When ah 2 > this line is the pre-image of the critical value. If this line intersects the critical circle, 
then the equilibrium is unstable and the line is the separatrix. If the line and the circle do not intersect, or 
if there is no circle (ah 2 < 0), then the line does not contain critical points, see Haller and Wiggins [2Uj . 
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4 Dynamics in two degrees of freedom 



To reconstruct the dynamics of H 3 , we have to attach a 1 -torus S 1 to every regular point of P£ . In this 
way the periodic orbits in P% give rise to invariant 2-tori and the unbounded reduced motions that exist for 
a = — 1 give rise to invariant cylinders; while regular equilibria lead to periodic orbits in K 4 . For the singular 
equilibria the isotropy group becomes important. The singular equilibria (tti, 7t 2 , 7r 3 , 7r 4 ) = (0, h 2 /(2a), 0, 0) G 
P£ 2 with /12 ^ lead to periodic orbits as well, albeit with half the period. From the singular equilibria 
7Ti = 0, we reconstruct the initial equilibrium at the origin of W 4 . 

The bifurcations of equilibria of the reduced system translate to bifurcations of periodic orbits in two 
degrees of freedom. The centre-saddle bifurcation turns straightforwardly into a periodic centre-saddle 
bifurcation. Where the singular equilibrium undergoes a Hamiltonian flip bifurcation we obtain a period- 
doubling bifurcation, supercritical (new orbit is elliptic) for a = +1 and subcritical (new orbit is hyperbolic) 
for a = — 1. Note that in the latter case the (hyperbolic) periodic orbit with twice the period bifurcates 
off from the Z 2 -isotropic one as hi increases through negative values (i.e. \h 2 \ decreases), while the elliptic 
periodic orbit with twice the period comes into existence as hi increases through positive values. 

In the latter case a = +1 of the 1:2 resonance the periodic orbits are elliptic and the regular fibres of the 
energy-momentum mapping EhA = (H 2 ,H 3 ) are compact. Whence the foliation into 2-tori is locally trivial 
and the pre-image £M _1 (U p) C M 4 of regular values p of £M is a 2-torus bundle. In the hyperbolic case the 
periodic orbits together with their (un)stable manifolds, determine the structure of this bundle. Together 
they turn the phase space K 4 into a ramified 2-torus bundle. A similar result for the case a = — 1 does not 
follow from general theorems. We therefore show by direct computation that the 1:— 2 resonance turns the 
phase space into a ramified cylinder bundle. With detuning there are some compact tori as well. 



4.1 Local triviality 

Because h — Ah + h 2 is linear in 7r 3 we can eliminate n 3 from the relation R} l2 = 0. To this end we define 
the polynomial 

Q(z) = 2(jc 2 z 2 (hi -z)- 4(A/i - \zf , (4.6) 

whence 7r| = Q(tv\) / '(4c 2 ). This is an elliptic curve, and can be parametrised by elliptic functions. Thus 
assume we have (iri(s; /), 7r 4 (s; /)) where s € R is the parameter and / = (h 2 , Ah) denotes the values of the 
energy-momentum mapping. By the foregoing analysis we know that the curve is not compact, but smooth 
when / is a regular value. By back-substiution we find 

7T2(S,/) = |(/l2-7Tl(s;/)) 

tt 3 (s; f) = ^(h-h 2 - Atti(s; /)) . 

Next we need to find a point (q,p) in phase space, so that the reduction mapping (q,p) i— > (tti, tti, 7T3, Wi) 
sends it to given values of the invariants. If such a mapping exists, we have a section of the bundle defined 
by the energy momentum mapping in a neighbourhood of a regular value. This mapping is 

I \ ( I — n ^3 7^4 

(Pi,qi,P2,qi) = V 71 "!; '— 

which is well defined because for regular values of / the invariant tti is nonzero (and positive in any case). 
Finally, since the S 1 acts on (p. q) by rotation, we have an explicit parametrisation of the cylinders in phase 
space given by 

(pi,qi,P2,q2)(s,4>;f) = 

(x/tti cos (j), yfni sin0, (7r 3 cos 2(f) + tt4 sin2</>)/7Ti, er(7r 3 sin 20 — tt 4 cos2</>)/7Ti) . 

The arguments (s; /) for all TTi have been suppressed. This is the local trivialisation of the bundle defined 
by the regular values of the energy momentum map. 
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4.2 Energy— momentum mapping 

The following results on the structure of the bifurcation diagram for the 1 : ±2 resonance were already 
described by Henrard [55]. However, his discussion of a = —1 was incomplete. Also see Schmidt and Sweet 
[5T] , Henrard [53J and Sweet [35] for a more general setting. 

Theorem 3. The critical values of the energy-momentum mapping are contained in the set given by the 
line h — h2 and the cubic curve 

[h 2 ,h-h 2 ] = |(f.-2A)[3t-2A,t 2 ]/c 2 . 

For t — the curve and the line h = h 2 are tangent at [h 2 ,h — h 2 ] = er[2A 2 , 0]/c 2 . 
For i = 4A/3 the curve has a cusp at [h 2 ,h- h 2 ] = -2cr[9A 2 , 8A 3 ]/27/c 2 . 
For t — 2A the curve intersects the line h — h 2 transversely in the origin. 



Proof. Eliminating tt^ in relation (3.51 by using the Hamiltonian gives the polynomial Q defined in the 
previous section in (4.6 1. Tangencies of the reduced Hamiltonian with the reduced phase space are given by 
the double roots of the cubic polynomial Q. Equating coefficients in —aQ(z)/(2c 2 ) — (z~r) 2 (z+at 2 /(2c 2 )) = 
gives the parametrisation of the cubic discriminant surface where ir% = r = at(t — 2A) /c 2 . A second solution 
branch with r = and h = h 2 for which h 2 — cr(4A 2 — t 2 )/2/c 2 gives a straight line of critical values. □ 

The characteristic feature of the curve of critical values of the energy-momentum mapping is that it 
always has a tangency with the straight line at t = 0. However, the cusp is below the line h = h 2 for aX > 0, 
and above for crA < 0. For A — the cusp is on the line. 

The transverse intersection at the origin marks the equilibrium point at the origin. The tangency is a 
bifurcation of periodic orbits that are relative equilibria of the reduced system. This bifurcation collides with 
the equilibrium point at the passage through resonance when A = 0. 

Because tti > the part of the discriminant locus given by t 6 [0, 2 A] is not part of the bifurcation 
diagram for a = 1; similarly for a = — 1 it is the complement. Since tt 2 = o(h 2 — 7Ti)/2 > and tti > only 
the part of h — h 2 with h 2 > is part of the bifurcation diagram; similarly for cr = — 1 it is the complement. 

For cr = 1 the piece of the straight line between the transverse intersection and the tangency is a stable 
relative equilibrium, beyond the tangency it is unstable. The type of the bifurcation is determined by the 
multiplicity of the preimage of the part of the curve beyond the tangency. The solution for iti when (h,g) 
are on the cubic curve is given by 

cr . „ XN (i-2A) 2 t-2A ,t-2\ 

? *(* " 2A )> — ^^ a ^' *i b -xr 

A point in the pre-image is given by pi = p 2 = q 2 = and 7Ti = q 2 /2. The flow of H 2 acting on this point 
gives all points in the pre-image. This is true for both signs of cr. Hence this bifurcation always is a period 
doubling bifurcation. 

The two branches emanating from the elliptic equilibrium (EE) at the origin are always elliptic (they 
are the nonlinear normal modes). For a — 1 one of them stays stable; while the other one loses stability in 
a supercritical period doubling bifurcation (supPD), and creates another elliptic relative equilibrium with 
twice the period. For cr = — 1 one of them disappears in the centre-saddle bifurcation (SC), while the other 
one collides with a hyperbolic relative equilibrium in a subcritical period doubling bifurcation (subPD) and 
only one hyperbolic relative equilibrium (with half the period of the previous hyperbolic relative equilibrium) 
remains. When A passes through zero, the role of the two elliptic branches emanating from the equilibrium 
is exchanged. 

Because a — — 1 gives a rotated critical curve the tangency now occurs for negative h 2 . The whole plane 
is in the image of the energy-momentum mapping when cr = — 1. When A = the only critical values are 
those on the ray h = h 2 with h 2 < 0. This is the case studied in |16[ 117] which has fractional monodromy. 
Our unfolding shows that at the endpoint of the ray "half a swallowtail" develops when the resonance is 
detuned. By continuity, encircling this whole structure and crossing the hyperbolic line h — h 2 for sufficiently 
large negative h 2 will give the same fractional monodromy. In conclusion we have shown 
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Figure 5: Image of the energy-momentum mapping for 1:2 resonance with c = 1, A = 1/2 (top left), A = (top 
right), A = —1/2 (bottom). Critical values are shown as thick red lines. Dotted lines are part of the discriminant 
locus, but are outside the image of the energy-momentum mapping (shaded), which is only to the right of the origin, 
bounded by the thick lines. The supercritical period doubling bifurcation is indicated by supPD. 
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Figure 6: Image of the energy-momentum mapping for the 1:— 2 resonance with c = 1, A = 1/2 (top left), A = (top 
right), A = — 1/2 (bottom). Critical values are shown as thick red lines. Dotted lines are part of the discriminant locus, 
but not critical values of the energy-momentum mapping. The whole plane is the image of the energy-momentum 
mapping. The preimage in the dark shaded region due to detuning lifts to one two-torus T 2 and a cylinder (figure 
|3] middle), whereas it is one cylinder otherwise (figure [3] top). Subcritical period doubling bifurcation, centre-saddle 
bifurcation, the elliptic equilibrium and the degenerate equilibrium are indicated by subPD, SC, EE and DE. 
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Theorem 4. The unfolding of the 1:2 resonance shows that nearby there is a supercritical period doubling 
bifurcation that passes through the equilibrium at resonance. The unfolding of the l:—2 resonance shows 
that nearby there is a subcritical period doubling bifurcation and a centre-saddle bifurcation that pass 
through the equilibrium at resonance. The unfolding of the l:—2 resonant equilibrium point has fractional 
monodromy. 

The proof of the last statement is deferred to Section [5j 
4.3 Perturbation analysis 

The equilibrium remains at the origin because the perturbation consists of higher order terms. Persistence 
of the periodic orbits and their bifurcations is straightforward. The Diophantine 2-tori will persist, provided 
that the Kolmogorov condition or iso-energetic non-degeneracy condition holds, see the next sections. For 
cr = —1 we will show that the iso-energetic non-degeneracy condition is violated along a line, while the 
Kolmogorov condition holds everywhere. Similar to the focus-focus case [TS], this can be considered to be 
an effect of the monodromy of the equilibrium. 



5 Monodromy 

The period of the reduced flow can be computed from 

TTl = {TVi,Tl(Tr)} = {iTl, Xni + C7T 3 } = 2c7T 4 . 



Fixing H = Ah and using (1.3 1, 7r| becomes a polynomial in tv\. Using dni — wdt the period integral is 
defined on the curve 

T = {(u>,7Ti) I w 2 = Q(vri) = 2cTC 2 TT 2 (h 2 - TTl) - 4(Ah - Xn^ 2 } . 

It is given by 

T{h 2l Ah) = j> 



w 



For the compact case a = 1 the meaning of this integral is clear. We now focus on the case a = — 1. A 
major problem in the non-compact case is to even define what is meant by period, rotation number, and 
action. Our approach is a simplified version of |29J . Assuming that there are no other critical points on the 
separatrix, we can split the dynamics into two parts: singular dynamics near the critical point and regular 
dynamics near the separatrix but away from the critical point, see figure [7] 
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Figure 7: Schematic representation of the dynamics for a — — 1 close to the separatrix of the singular point at the 
origin. The green dashed curve indicates the continuation of the stable and unstable manifold of the singular point 
when compactified. 

Assuming that the system has compact invariant tori near the separatrix the contribution of the dynamics 
far away from the equilibrium is some bounded function with smooth dependence on the initial conditions 
as long as we stay on the same "side" of the separatrix. This function contains the symplectic invariants 
of [25], but we are going to completely neglect it. As in 29J, the closer we get to the critical value, the 
main contribution to the period comes from dynamics near the critical point. This contribution diverges 
in a characteristic way. It is this universal behaviour that we are interested in. For a hyperbolic point the 
period diverges logarithmically. Our equilibrium is not hyperbolic. Instead we find that the period diverges 
algebraically. There is one additional trick, which has first been used in [T2J . Instead of integrating along the 
orbit until it leaves a neighbourhood of the critical point, we integrate until the orbit reaches infinity. Thus 
the integral involved becomes nice namely, a complete elliptic integral. The difference between integrating 
up to some arbitrary cut-off point and integrating to infinity is a bounded smooth function, which can be 
neglected. This becomes clear when considered from the point of view of algebraic geometry. The frequency 
map is given by integrals of holomorphic one-forms over cycles on a complex affinc cllpitic curve. After 
adding the point at infinity to make the affine elliptic curve a compact elliptic curve in projective space, 
cycle that go to infinity on the affine curve are now compact. After this compactification the integral of 
the holomorphic one form over the compact cycle gives the frequency map. The contribution near the point 
at infinity is smooth. The asymptotics of the frequency map near the equilibrium point is obtained from 
the contirbution of the cycle near the equilibrium point. This behaviour is unchanged by adding higher 
order terms to the Hamiltonian. This is the reason why our analysis gives the asymptotic behaviour for all 
completely integrable systems with the same normal form up to order three. 

The parameter c can be removed by the scaling (tti, h 2 , Ah) — > 1/c (iri, h 2 , Ah) (similar to the scaling 
in section 2). The cycle 7 extends around the largest real root of w and infinity. This root has the same 
sign as Ah. For values (h 2 , Ah) near the line h — h 2 there are three real roots. The period integral is finite 
everywhere except when approaching the critical values along the line h — h 2 , T(/i 2 ,0) = 00 for h 2 < 0. 
Even though the reduced motion is unbounded, the time it takes to reach 00 is finite. This explains why we 
can treat everything formally as if we were in the compact case. 

To compute the rotation number we can cither differentiate the action A (sec below) with respect to h 2 , 
or derive the differential equation for dynamics in the fibre. We will pursue the latter. 

The 5' 1 -action makes each cylinder a principal <5 -bundle, with vertical direction along the fibre. Singling 
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out a horizontal direction amounts to defining a connection. While an angle in the fibre given by 



= arg z\ 

is not well-defined, the 1-form d0 does define a principal connection simultaneously on all cylinders. The 
formula {0, H 2 } = 1, however, can still can be used since it only involves derivatives. The Poisson brackets 
with 0, which describe the dynamics in the fibre in phase space are computed to be 

{0,r/} = 1, {0,^} = 1, {0,7r 3 } = 7r 3 /7ri. 

Thus the dynamics of generated by TL is given by 

„, r„ , „ , 7Tq Ah — A7Ti Ah 
0' = {Q,H} = {0, Atti + C7r 3 } = A+c— = A + -= — . 

7Tl 7Tl 7Ti 

For later use we define an integral that describes the varation of when this dynamics in the base completes 
a full cycle. Changing the integration from dt to d7ri and integrating over a period T gives the winding 
number as 

W(h 2 ,Ah) = ^I^^. (5.7) 

27T J 7 7Ti 10 

We call this integral the winding number, because it measures the (fractional) monodromy of the systcmf^] 
Note that this is not the rotation number of 7i 3 , because H. and 7i 3 differ by H 2 . This winding number 



integral ( 5.7 ) has a pole at 7Tl = 0. The residue of this pole is j^i sgn Ah. This integral shares many properties 
with the one discussed in |17j . However, our integral only describes the dynamics near the resonance, while 
the integral in |17j the terms added for compactification also play a role. The singular behaviour near the 
equilibrium, however, is the same. 



As already mentioned, for the computation of the rotation number of the full Hamiltonian (3.4 1 the term 
H 2 needs to be included. The rotation number of the Hamiltonian without this term is what we defined to 
be the winding number in the previous paragraph. The equations of motion of the two Hamiltonians differ 
by 1, which originates from {0,_ff2} = 1- Therefore the dynamics of generated by the full Hamiltonian 
H3 is given by 

= {0,H 3 } = 1 + — . 

Changing the integration from dt to d7ri and integrating over a period T then gives the rotation number 

2nR(h 2 ,Ah) = T(h 2 ,Ah) + 2nW(h 2 , Ah) . 

The period integral has logarithmic singularities. Therefore R does as well. For W the situation is not so 
obvious. 

In order to understand the behaviour of the integrals near the origin we perform the scaling iri = g 2 z 
and introduce scaled parameters 

(h 2 ,Ah, A) -> (g 2 h 2l g 3 Ah, gX) . 

This is similar to the scaling done in section 2, but we repeat these steps here because we want to interpret 
this scaling as choosing a new co-ordinate system on the image of the energy-momentum map. This scaling 
shows that the winding number is independent of the "radius" g. In order to study the behaviour of W when 
the origin is encircled we therefore introduce weighted polar co-ordinates in the image of the momentum 
map by 

(h 2 ,Ah,X) = (g 2 cos9,g 3 sm9,gX) . 

Even though a winding number is usually an integer, we will use the term for the function W that changes by a (half-) 
integer over a cycle around the origin of the energy-momentum map. 
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The transformation from (/12, Ah) to new variables near the origin of the momentum map (u, v) = £> 5/,2 (cos 9, sin 9) 
is C° at the origin and C°° everywhere else. With this notation the crucial integrals become 

T(h 2 ,Ah) = -A(9), 2nW(h 2 ,Ah) = B(0). (5.8) 
Q 

If the degenerate equilibrium is approached on a curve in parameter space {Ji2,Ah) with non- vanishing 
derivative at the origin, the leading order term of T reads | Ah\ 1 . Along the line Ah = 0, h 2 > 0, we find 
l/\fh~2~ as leading order term of T. So in our two-degree of freedom system the leading order behaviour of 
the period T as a function of (h 2 , Ah) depends on how the equilibrium is approached. This is in contrast to 
one degree of freedom systems where the energy h is the only parameter. 
The integrals A and B are defined on the curve 

f = {0, w) : w 2 = 2az 2 (cos 9 - z) - 4(sin 9 - \z) 2 } 

where 

dA(9) = \dz, dB{9) = ^dz . 
w zw 

The parameter A must be chosen sufficiently small so that there are no critical values encountered on the 
loop around the origin in the (h 2 , Ah) plane. The plots of the functions A and B (i.e. of T and W on a loop 
around the origin with g — const) are shown in figure [8] 

The period T diverges algebraically as \ j q when the origin is approached. This is due to the fact that 
the equilibrium is a degenerate singularity: specifically, the value of the Hessian of the cubic integral at 
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the critical point (i.e. the origin) vanishes identially. The winding number is independent of the distance 
from the origin g. The period diverges logarithmically when 9 — > ±tt. This is the usual divergence when 
approaching a hyperbolic periodic orbit. The winding number is only a continuous function of 9. It is not 
differentiable. Surprisingly, the winding number is not a periodic function of 9, even though its defining 
integrand is. 

We now prove that the increase in values of W (or B) when we let 9 G [— 7r, tt] is indeed as sugested by 
figure [8] First of all observe that B is an odd function, so we can restrict to 9 > 0. For A = the polynomial 
w 2 has one non-negative real root and two roots with non-positive real part (which may be real or complex) . 
For 9 = these are (1, 0, 0) and for 9 = tt these are (0, 0, —1). There is a double root for 9 = 0. But since we 
always integrate from the largest root onwards this does not concern us (there is another collision of roots, 
but again it is not in the physically allowed region). Thus the function B(9) is smooth for 9 E (— tt, tt). The 
double root at for 9 = tt is at the boundary of the integration interval, so care needs to be taken when 
evaluating B. At first sight it looks like B(tt) = (which would also be required for a periodic odd function). 
However, we will now show that B{9) is not periodic, but instead has non-zero tt/2 value at 9 = tt. Since 
B is odd therefore B(—ir) = —tt/2 and B(tt) — B(—tt) = tt: the value of B increases by tt when 9 traverses 

[-7T, TT]. 

Lemma 5. 

lim B(9) = ±tt/2 + arcsin (V2\) 
e^±7r=F V / 

Proof. We consider the limit to +tt from the left first. When 9 — > tt we can treat 5 = sin 9 as a small 
parameter. The colliding roots behave like <5, thus we scale z = yS and find 

dB = l/(y*J-2y 2 (cos9- y$in9) - 4(1 - \y) 2 )dy 

The root that limits to z = — 1 in the new variables tends to -co when 5 — * 0. The other two roots of 
opposite sign tend to non-zero values. Denote these roots (for any small value of S) by yoo < y~ < y+. The 
integration interval now is [y + ,oo). In order to regularize the square root singularity of the integrand at y + 
and also make the lower integration boundary independent of S we introduce new variables by y = y+ + u 2 . 
Now the integrand is a smooth function of u and 5 uniformly over the integration interval. We now can 
compute the Taylor series in 9 — tt. The first term gives the value of B at 9 = tt as 

V2du 

(u 2 + y+)^/2y + + 4A(1 - Ay+) + u 2 (l - Wj 

Since the integrand is even function of u, the interval of integration can be extended to [-co, oo], absorbing 
one factor of 2. Direct evaluation of the integral gives the stated result. The computation is similar for 
B(tt), except that now d is negative which changes the sign of the constant term in the result, but not the 
A-dependent part. As a result the difference B{tt) — B(—tt) = tt is independent of A. □ 

For later use we also record the following 

Lemma 6. The function B(9) for A = is strictly monotonic. 

Discussion In the appendix we give the expression for B' in terms of Legendres complete standard integrals 
K and E. Using these formulas it can be numerically unambiguously verified that B' does not have a zero. 
The function B' can be plotted with any standard symbolic algebra program. We find B' > 0.15, so it is 
bounded away from zero such that only a few digits of accuracy are needed. Thus round-off errors in the 
computation do not play a role. We have not been able to find an analytic proof yet. end of discussion 
To establish monodromy is now simply a matter of finding an expression for the actions / = (/i,^) for 
the integrable system. One action I2 is trivial, namely the generator of the S 1 symmetry. Set I2 = /12 ■ The 
non-trivial action I\ of the system (or more precisely, its leading order term near the equilibrium) contains 
some terms that depend smoothly and periodically on 9. However, it has a term Wh 2 , which changes by 
1 over a cycle around the origin, see Lemma |E] Thus 1\ changes by 1 12 after traversing this loop. This 
behaviour is called fractional monodromy [27 . 



B(tt) = 2 / 
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Lemma 7. The non-trivial action I\ of the l:—2 resonance has the form 



n A I 

I 1 (ha,Ah) = -—T-k 2 W. 

I Z7T 

It is discontinuous across the singular value at Ah — and it changes by hill — I2/2 along a loop around 
the origin. 

Proof. For the computation of the action of the reduced system two canonically conjugate functions whose 
Poisson bracket equals one are needed. It is easy to check that 

{COS —==,7r 2 } = l 
V7rf7T 2 

before introducing 77 in the Poisson structure. Now the action is given by 

Integration by parts and changing the integration variable to ni = h-z + 2iT2 gives 
Ah f f 3 h 2 \ diti 



h(h 2 ,Ah) 

ZlT J \Z 7Tl J w 

and the result follows. □ 

The formula of Lemma [7] should be compared to the analogous expression for a focus- focus point, which 
is I\ = J\T — J 2 W where Ji are the co-ordinates from the Eliasson normal form [29 . The two formulas are 
surprisingly similar, even though the corresponding equilibria are quite different. 

Note that the remarkable formula of Lemma [7] can also be read as a decomposition of the angle change 
2~kW into a geometric phase (namely the action) and a dynamic phase (the term proportional to T) . In many 
cases simple relations like this can be found for the change of angle in the fibre computed from quantities in 
the base, see, e.g., Montgomery [25] . 

For us the most important aspect of Lemma [7] is that together with Lemma [5] it immediately implies the 
following monodromy theorem, complementing the work of [271 117j . 

Theorem 8. The Hamiltonian system describing the l:—2 resonance has fractional monodromy |. 

Let us point out the philosophy in our approach once more. Following Duistermaat [10] classical mon- 
odromy can be understood by considering the period lattice, or cquivalently, following Cushman [9] the 
rotation number of the system. Following Zhilinskif [35j another way to think of Hamiltonian monodromy 
is in terms of lattice defects in the associated quantum analog. In the present case the rotation number 
is singular at the critical values of the energy- momentum map along the line h = h 2 , and the fibres are 
unbounded. This makes it hard to speak about continuing the period lattice or rotation number around the 
origin (the degenerate critical point corresponding to the equilibrium) . The first papers discovering fractional 
monodromy [27] used a geometric approach. Later [16j a rescaled time was used in order to "remove" the 
blowup of T. Our approach is similar to this, in that we are dealing with a similar type of elliptic integrals. 
There are two differences. Firstly, in the proof of monodromy we directly consider the actions of the system. 
Even when the rotation number is singular, the action can always be chosen continuous across separatrices 
(up to integer factors counting the number of tori, etc.). In our case the action can be made continuous at 
any point, but globally it is not single valued, due to the contribution of W . It seems that in general the 
actions are the more "robust" objects that are well suited for this type of consideration. Secondly, instead 
of compactifying the Hamiltonian by adding higher order terms we directly treat the non-compact situation, 
by only considering dynamics near the singularity. This is inspired by the approach of San Vu Ngoc [29] to 
the symplectic invariants near focus-focus points. Our main new results are presented in the next section, 
where we consider the derivatives of the functions A and B in order to study the non-degeneracy condition 
of the KAM theorem. 
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Figure 9: Illustration of the condition on vanishing twist in the 1 : —2 resonance. The twist condition is violated 
where the lines of constant h are tangent to the contour lines of R(fi2, Ah) which is the case for Ah — and h^ > 0. 
h-2 = Q 3 cos 9. 

6 Frequency Map 

The results in this and the following section are valid at the bifurcation A = only. For perturbation theory 
we have to consider the actual rotation number R not just the less singular winding number part. To study 
the derivatives of R and the frequencies it is best to use the scaled forms of the integrals. Only the leading 
term when g — > is relevant. Inverting d(h 2 , h)/d(g, 9) gives 

= -Ag 2 sm9 + 0(1/ q), ^- = -2Agcos 9 + 0(1/ g 2 ) , 
on 2 Cft-2 

where the determinant of the transformation is 

2 

£> 4 (5-cos26») ' 

Theorem 9. ('Vanishing Twist) The isoenergetic non-degeneracy condition is violated on a line emerging 
from the critical value at the origin of the l:—2 resonance. Namely the function R(h^) has a critical value 
at h 2 = h for A = 0. 

Proof. Computing the derivative of 2irR — A(9)/g+ B(9) after scaling gives 

27r !^ = -4#+ (~ + B ') |^ = A(A sin 0-2,4' cos 0) + O(l/g 3 ) 
dh 2 g 2 oh 2 \ g J oh 2 

For A = the function A is even in 9 and therefore the leading order term vanishes at 9 — 0. For A ^ the 
origin becomes a non-degenerate elliptic fixed point, and the vanishing twist will instead be present near the 
centre-saddle bifurcation [T2] , see the bifurcation diagram figure [6j □ 

A graphical illustration of the condition on vanishing twist is given in figure [9] 
The frequencies of the system are 

uii = 2vr/T, and uj 2 = 2irR/T = 1 + 2itW/T . 



d(g,9) 
d(h 2 ,h) 
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Figure 10: The non-violation of the Kolmogorov condition means the frequencies plotted in the image of the 
energy-momentum mapping give a "nice" co-ordinate system. The contour lines of u>2 are black, the ones of u>i gray. 



Theorem 10. The Kolmogorov non-degeneracy condition 
dl T 

is satisfied at every point in a neighbourhood of the origin of the l:—2 resonance. 
Proof. The determinant is the product of three determinants 



dl 



d(Q,9) 



d(h 2 ,h) 



d{h 2 ,h) 



The middle one is A, and hence non-zero outside the origin. Since h 2 — I 2 the last one reduces to —dh/dli = 
—oji which is non-zero outside the critical values. The non-trivial first determinant is (exactly!) given by 



doj 




d( B ,9) 


= ( 



-2itqA'/A 2 



2lTQ 



B' 
A 2 



In Lemma [6] we have shown that B(9) is monotonic, see figure [8j This completes the proof. 
Collecting the above results we find that 

Notice that this diverges like 1/q 2 when approaching the origin away from 8 — ±7r. When 9 — > ±7r for fixed 
p, however, the determinant approaches zero because both, A and B' diverge logarithmically. □ 

In this section we only considered the case without detuning, A = 0. However, using the analysis for 
period doubling in [21] and for the centre-saddle bifurcation in [12] it is straightforward to get a qualitative 
understanding of vanishing twist for A ^ as well. For a = +1 no vanishing twist is forced by the bifurcation. 
For a = — 1 the cusp of the centre-saddle bifurcation will produce a curve emanating from it along which 
the twist condition is violated. For A — > this curve will locally coincide with A/i = as predicted here. 
Combining the images of the energy-momentum map for different parameters A into a three-dimensional 
picture there will therefore be a surface along which the twist condition is violated. To our knowledge the 
Kolmogorov condition has not been investigated near these bifurcations. 
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A Standard Form of the Integrals A, B, A', and B' 

Here we present the main elliptic integrals in Lcgcndrc normal form for the cases of three real and one 
real root of w 2 . This can be avoided using a Landen transformation [5] but since this doesn't simplify the 
Lcgcndrc normal forms we refrained from doing so. We denote the real roots by a > b > c and the complex 
conjugate ones by b and b. 
A(9) is simply given by 



A(9) = V2gK{ 7 



(A.9) 



where 



9 



yja — c a — c 

for three real roots and 

1 d + bi — a 

Vd 2d 



d 2 = (Rc6-a) 2 + (Im&) 5 



for one real root. The derivative of A is of 2nd kind. For three real roots we find 
A' (6) = V2g [K(m) Oi + au 2 ) - u 2 (a - c)E(m)} 

with 

sin 29 



(A.10) 



(/x 2 cos 



0-9), u 2 



jJL 2 COt ( 



(cos 2(9 -5) , 



2D ^ " AD 

where D is the discriminant of uj 2 given by D — 2/j 2 cos 3 9 + 27 sin 2 9. In the case of one real root we have 
A' (9) = V2g[K(m) ( Ul + (a + d)u 2 ) -2du 2 E(m)} . (A.ll) 
B{9) is the elliptic integral of 3rd kind 



B(9) = V2g sin (9 
e case of thr 
B{9) = V2- 



K{m) 



b ba 
in the case of three real roots and 

K(m) 



n ( — , m 
a 



a(l — a 2 ) 



n 



(A.12) 



(A. 13) 



otherwise, where 
a — d 

a = 



a + d 

Its derivative is of 2nd kind. For three real roots we find 
B'(9) = -V2g [K(m) (m + au 2 ) - u 2 (a - c)E(m)} 

where 

u 2 cos9 . „. 3/i 2 , . o „\ 

ui = f -^-(co S 29-5), u 2 = ^(2 + sm 2 9). 

In the case of one real root we finally have 

B'(9) = -V2g [K{m) (u x + (a + d) u 2 ) - 2 d u 2 E[m)] 



(A.14) 



(A. 15) 
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